Preprint typeset using I4TgX style emulateapj v. 08/13/06 



IS GRAVITATIONAL LENSING BY INTERCLUSTER FILAMENTS ALWAYS NEGLIGIBLE? 

Martin Feix', Dong Xu^, HuanYuan Shan\ Benoit Famaey*, Marceau Limousin", HongSheng Zhao'-^ and Andy Taylor^ 

ABSTRACT 

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), jn ~ 10 - 
10"^. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, 
a single filament inclined by x 45° from the line of sight can cause substantial distortion of background sources 
pointing towards the filament's axis {k = j = {1 - A"')/2 ~ 0.01); this is rigorous for infinitely long uniform 
filaments, but also qualitatively true for short filaments (~ 30Mpc), and even in regions where the projected 
matter density of the filament is equal to zero. Since galaxies and galaxy clusters are generally embedded 
in filaments or are projected on such structures, this contribution complicates the interpretation of the weak 
lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing 
order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous 
weak lensing signals, e.g. the "bullet cluster" of Clowe et al., the "cosmic train wreck" of Abell 520 from 
Mahdavi et al., and the "dark clusters" of Erben et al., filamentary structures might contribute in a significant 
and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing 
signal could, in principle, be used to falsify MOND/TeVeS and its variations. 

Subject headings: cosmology: theory - dark matter - gravitation - gravitational lensing - large-scale structure 
of the universe 



1. EVITRODUCTION 

Without resorting to cold dark matter (CDM), the modi- 



fied Newtonian dynamics (MOND) paradigm ( Milgrom 1983 



Bekenstein & Milgrom 1984) is known to reproduce galaxy 
scaling relations like the Tully- Fisher relation (|Tully & Fisher] 
|1977| l, the Faber-Jacks on law (|Faber & Jackson||1976[ ) and 
the fundamental plane ( |Djorgovski & Davis|1987| l) as well as 
the ro tation curv es of individual galaxies over five decades in 
mass (Sanders & McGaugh 2002, Bekenstein 2006; Sanders 



& Noo rdermeer 2007; Milgrom & S anders|2007| |2003 [ |Nipoti| 
et al. ]2007, Famaey et al. 2007 af! In particular, the re- 



cent kinematic analysis of tidal dwarf galaxies by |Boumaud| 



et al. (2007 1 is very hard to explain within the classical CDM 
framework while it is in accordance with MOND ( |Gentile| 



et al. 2007; Milgrom 2007) 1. In addition, observations of a 
tight correlation between the mass profiles of baryonic mat- 
ter and dark matter in relatively isolated (field) galaxies at all 
radii ( |McGaugh et al. 2007 Famaey et al. 2007b j ) are most 
often interpreted as supporting MOND. Nevertheless, in rich 
clusters of galaxies, the MOND prescription is not enough 
to explain the observed discrepancy between visible and dy- 



namical mass (jSanders 2007; Angus et al. 2007; Feix et al. 
2008| l. At very large radii, the discrepancy is about a factor of 



2, meaning that there should be as much dark matter (mainly 
in the central parts) as observed baryons in MOND clusters. 
One solution is that neutrinos have a mass at the limit of detec- 
tion, i.e. ~ 2 eV, which can solve the bulk of the problem of 
the missing mass in galaxy clusters, but other issues remain 
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(Angus et al. 2008] l. These 2 eV neutrinos have also been 
invoked to fit the angular power spectrum of the cosmic mi- 
crowave background (CMB) in relativistic MOND (Skordis] 
|et al. 2006), and are thus part of the only consistent MOND 
cosmology presented so far. In the following, we will refer to 
this model as the MOND hot dark matter (yuHDM) cosmology 
( [Angus et al.|2007j i''. 

On the other hand, strange features have recently been dis- 
covered in galaxy clusters, which are hard to explain, such as 
the "dark matter core" devoid of galaxies at the c enter of the 
"cosmic train wreck" cluster A bell 520 (M ahdavi et al.|2007| l 
and others ( |Erben et aL]|2000[ [Clowe et al.||2006i ). Here, we 
consider the possibiUty that this kind of features could be due 
to the gravitational lensing effects generated by an interclus- 
ter filament in a universe based on tensor- vector-scalar gravity 
(TeVeS; Bekenstein 2004 ), one possible relativistic extension 
of MOND (cf. 'Zlosn ik et aL|2006| [2007) |Zhao|2007| ). How- 
ever, we are not performing a detailed lensing analysis of any 
particular cluster in the presence of filaments, but rather pro- 
vide a proof of concept that the influence of filaments could 
be much less negligible in a MONDian universe than within 
the framework of general relativity (GR). 

Filaments are among the most prominent large-scale struc- 
tures of the universe. From simulations in ACDM cosmology, 
we know that almost every two neighboring clusters are con- 
nected by a straight filament with a length of approximately 
20 - 30 Mpc ( [Colberg et al.|2005 ). For instance, the dynamics 
of field galaxies, which are generally embedded in such fila- 
ments, as well as their weak lensing properties are persistently 
influenced by this kind of structures, generally encountering 
accelerations of about 0.01 -0.1 x 10" m s"^. Filaments also 
cover a fair fraction of the sky, much larger than the covering 
factor of galaxy clusters. Thus, there is a good chance that fil- 

* Note, howev er, that one could also switch to steri le neutrin os with masses 
of a fe w eV (e.g. |Gentile et al.|2008|rMaltoni et al.] 2004 Malto ni & Schwetz| 
[2007^ and that massive (sterile) neutrinos are not indispensable within cer- 
tain covariant formulations of modified gravity, e.g. the I'A model, which 
can mimic the effects of neutrinos in clusters and cosmology as well as the 
behavior of a cosmological constant jZhao|2007| . 
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aments might be superimposed with other objects on a given 
Une of sight, hence affecting the analysis of observational data 
like, for example, weak lensing shear measurements. Such re- 
cent studies prompted us to investigate the possibility that, in 
the context of MOND, end-on filamentary structures could be 
responsible for creating anomalous features in reconstructions 
of weak lensing convergence maps such as the peculiar "dark 
matter core" devoid of galaxies in Abell 520 (Mahdavi et al. 
[2007] l. 

Short straight filaments are structures which, at the best, 
are partially virialized in two directions perpendicular to their 
axis. According to |Colberg et al.| ( |2005[ ), a filament gener- 
ally corresponds to an overdensity of about 10 - 30, having a 
cigar-like shape. Furthermore, filamentary structures tend to 
have a low-density gradient along their axis and, in the per- 
pendicular directions, they have a nearly uniform core which 
tapers to zero at larger radii, usually about 2-5 times their 
core radius. Since filaments are typically much longer than 
their diameter, we shall approximately treat them as infinite 
uniform cylinders of radius Rf - 2.5 hr^ Mpc. 

Lacking a MOND/TeVeS structure formation A^-body sim- 
ulation (with or without substantially massive neutrinos), we 
shall adopt the naive assumption that filamentary structures 
have roughly the same properties in MOND and in CDM, 
which will be justified in ^ Deriving expressions for the 
TeVeS deflection angle and setting up a cosmological back- 
ground, we conclude that the order of magnitude of the TeVeS 
lensing signal caused by filaments is compatible with that of 
the previously mentioned observed anomalous systems. In 
addition, we find that there is fundamental difference between 
GR and MOND/TeVeS for cylindrically symmetric lens ge- 
ometries (see Fig. [TJ; in contrast to GR, the framework of 
MOND/TeVeS allows us to have image distortion and am- 
plification effects where the projected matter density is equal 
to zero. As for a more realistic approach, we also consider 
a model where the filament has a fluctuating density profile 
perpendicular to its axis. Compared to the uniform model, we 
find that the lensing signal in this case is smaller, but still of 
the same order, taking into account that the filamentary struc- 
tures may be inclined to the line of sight by rather small angles 
{9 < 20°). Finally, we demonstrate the impact of filaments 
onto the convergence map of other objects by considering su- 
perposition of such structures with a toy cluster along the line 
of sight. Again, our results show an additional contribution 
comparable to that of a single isolated filament. 

2. MODELING A FILAMENTARY LENS 

We investigate the effect of gravitational lensing caused by 
a straight filament connecting two galaxy clusters in both GR 
and TeVeS gravity, henceforth using units with c - I. As 
a first simple approach, we shall take the filament's matter 
density profile to equal an infinitely elongated and uniform 
cylinder which is illustrated in Figure [T] The cylinder's line 
density. 



A = MIL = pnR 



(1) 



is taken to be constant, where M is the total mass, L denotes 
the length along the symmetry axis, /?/ is the cylinder's radius, 
and p is the volume density. A photon traveling perpendicu- 
lar to the filament's axis will change its propagation direction 
when passing by the cylinder due to the local gravitational 
field which is assumed to be a weak perturbation to flat space- 
time, i.e. all further calculations may be carried out within the 
non-relativistic approximation. In this case, it is well-known 
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Fig. 1. — Light deflection by an infinitely elongated cylinder of constant 
mass density; the unperturbed photon traveling along the z-direction passes 
the filament at the distance y (impact parameter) from the filament's axis and 
is deflected by the angle a. The line density of the filament is assumed to 
be constant, A = M/L = pnR^, where p is the volume density and Rf is the 
cylinder's radius. 

( Schneider et al.|1992 i that the photon's deflection angle can 
be expressed as 



(2) 



where €) is the total gravitational potential, Vj^ denotes 
the two-dimensional gradient operator perpendicular to light 
propagation and integration is performed along the unper- 
turbed light path (Born's approximation). In our example (see 
Fig. [T]), the filament's axis is aligned with the jc-axis, and light 
rays propagating along the z-direction are dragged into the in- 
directions due to the symmetry of the resulting gravitational 
field. Keeping this configuration and introducing cylindrical 
coordinates, we may rewrite equation (|2| as 



oo 

r o' 

a(y) = 4y dr. 



(3) 



where the prime denotes the derivative with respect to the 
cylindrical radial coordinate r, i.e. A - dA/dr. Note that 
even in the context of MOND/TeVeS, we may still assume that 
most of the light bending occurs within a small range around 
the lens compared to the distances between lens and source 
and observer and source, thus enabling us to fully adopt the 
GR lensing formalism. 

In gravitational lensi ng, it is convenient to i ntroduce the de- 
flection potential *P(0) ( Schneider et al. 1992 1: 



*I'(^) = 2 



r o(A^, 



z)dz. 



(4) 



where we have used 6 - ^/Di. Here ^ is the two-dimensional 
position vector in the lens plane, and Ds, Di, and D/, are 
the (angular diameter) distances between source and observer, 
lens and observer, and lens and source, respectively. If a 
source is much smaller than the angular scale on which the 
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lens properties change, the lens mapping can locally be lin- 
earized. Thus, the distortion of an image can be described by 
the Jacobian matrix 



1-K-yi -j2 
-72 I- K + Jl 



(5) 



where P=f]/Ds and if denotes the 2-dimensional position of 
the source. The convergence k is directly related to the deflec- 
tion potential through 



1 



and the shear components yi and 72 are given by 



1 



.del 



dol, 



72 = 



801802 



(6) 



(7) 



Because of Liouville's theorem, gravitational lensing pre- 
serves the surface brightness, but it changes the apparent solid 
angle of a source. The resulting flux ratio between image and 
source can be expressed in terms of the amplification A, 



(8) 



Considering the symmetry properties of our cylindrical lens 
model and the configuration in Figure [T| equation (|6| further 
simplifies to 



^ ^ 1 D,Di, da(y) 



(9) 



with the convergence k being related to the quantities y (7^ - 



j\ and 72 = 0) and A as follows: 



1 -A- 



K — y — 



(10) 



Furthermore, let us introduce the complex reduced shear g 
given by 

n + '72 



l-K 



(11) 



This quantity is the expectation value of the ellipticity x of 
galaxies weakly distorted by the lensing effect, thus corre- 
sponding to the signal which can actually be observed. In 
our case, we find that the absolute value of the reduced shear 
is 1^1 - 7/(1 - k), and assuming that k - y <& 1, we obtain 
\g\ ~ K - 7. Note that the above result is independent of the 
particular law of gravity. 



Inserting the above into equation (|9]), we may obtain the cor- 
responding convergence field. As expected, kn equals zero 
outside the cylinder's projected matter density. 

(II) For y < Rf, the deflection angle has to be calculated 
from 
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(14) 



Carrying out the integrations in equation ( [T4] i, we finally end 
up with the following expression: 



ONiy) = 



2GA 



^ . ■ y 
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Rl \R 



f 



(15) 



Using equation (|9]l, the convergence in this case turns out to 
be 

'<N(y)-2^^M^^. (16) 

nRj ' 

2.2. MONDian Case 

Now we shall consider light deflection within the frame- 
work of TeVeS gravity. Choosing a certain smooth form of 
the free interpolating function yu (fo r further details see,Beken- 
stein 2004) which has been used by Zhao et al. ( 2006| l, the to- 
tal gravitational acceleration may be written in the following 
way: 



8M(r) = |VOM(r)| = g^ir) + y/gN(r)ao, (17) 



with r again being the cylindrical radial coordinate and ^M(r) 
the total non-relativistic gravitational potential in TeVeS. The 
constant ao = 1.2 x 10"'" m s"^ characterizes the accelera- 
tion scale at which MONDian effects start to become impor- 
tant compared to Newtonian contributions. Since filaments 
are the most low-density structures within the universe, their 
internal (Newtonian) gravity is very small. Therefore, the 
MONDian influence yields an enhancement of the gravita- 
tional field which is on the order of ao/gN, being extremely 
large in such objects. For this reason, we may expect a sub- 
stantial difference concerning the lensing signal caused by fil- 
amentary structures in TeVeS. Equipped with equations ([3]), 
( [T2] i and ( [T7| ) we are ready to proceed with the analysis of our 
cylindrical filament model: 
(I) For Rf <y, the deflection angle is given by 



2.1. Newtonian Case 

The Newtonian gravitational field of our filament model is 
given by 



gNir) = \V<S>N(r)\ = 



GA r 
2nR 



2' '•<^/' 

/ 



(12) 



GA 1 
2n r ' 



r>R 



with A being the previously defined line density given by 
equation ([T]i. 
(I) For Rf <y, evaluating integral (|3]l yields 



ctN(y^ - GA - const. 



(13) 



^^GAao^ J 



dr 



^J2GAaoy. 



r(3/4) 

In this case, the convergence reads as follows: 



KM(y) 



D/D/,r(l/4) GAap 
D, r(3/4)V Sy 



(18) 



(19) 



(II) For y < Rf, integral (|3]l has to be split in several parts, 
similarly to equation ([T4|. Using elementary calculus, we fi- 
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Fig. 2. — Density profile p(r) {top left), radial evolution g(r) of the total gravitational acceleration {top right), deflection angle &{y) {bottom left) and convergence 
K{y) {bottom right; a: = y = (1 - A"')/2) in Newtonian {dashed line) and IVIONDian (solid line) gravity for th e un iform filament cylinder model whose axis is 
inclined by an angle 6 = 90° to the line of sight, assuming zi = 1, Zs = 3, and the flat pHDM cosmology of eq. j23j in IVIOND/TeVeS. The radius of the filament 
is Rf = 2.5 /i"' IVIpc, and the overdensity within the filament is taken as 20 times the mean density po in accordance with the results of Colberg et al. (2005). 
Note that, for consistency, the Newtonian results are based on a flat ACDM cosmology with 0.,„ = 0.3 and = 0.7. 



nally arrive at 



3 1 
4' 2 



(20) 
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where UNiy) is the Newtonian deflection angle given by equa- 
tion ( [TS] ) and 3(p q)(a,b) is the generalized incomplete Beta 
function defined by 

S(p,g-)(a,b) = J f'-\l -tf-^dt, Re(a),Re(b)>0. (21) 
p 

As the expression for the convergence km turns out to be quite 
lengthy, we will drop it at this point. 

From equations ([T8| and ([19]), we find that the deflection 
angle aM outside the cylinder's projection increases with the 



square root of the impact parameter y {un - const), the 
convergence km decreases with the inverse square root of y 
(Kff - 0). In fact, this reveals quite a fundamental difference 
between MOND/TeVeS and GR; since - 0, we also have 
jN = and An - I according to equation ( [TO| i, meaning 
that there will be no distortion effects, as well as no change 
in the total flux between source and image, i.e. wherever the 
projected matter density is zero, the lens mapping will turn 
into identity. However, this is no longer true in the context of 
MOND/TeVeS as the convergence and the shear field do not 
vanish (see Fig. |2]). Obviously, this is a case where the MON- 
Dian influence does not only enhance eff'ects that are already 
present in GR, but rather creates something new, which, in 
principle, could be used to distinguish between laws of grav- 
ity (see 

Varying the inclination angle 6 of the filament's axis to the 
line of sight, the lensing properties derived in this section have 
to be rescaled by a factor of sin ' 9 in both GR and TeVeS. 

3. MODEL APPLICATION 
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Fig. 3. — Density profile p(r) {top left), radial evolution g(r) of the total gravitational acceleration {top right), deflection angle &(>■) {bottom left) and convergence 
Kiy) {bottom right; a: = y = (1 - A"')/2) in Newtonian {dashed line) and MONDian (solid line) gravity for the unifor m fil ament cylinder model whose axis is 
inclined by an angle d = 90° to the line of sight, assuming z; = 1, Zj = 3, and the flat minimal-matter cosmology of eq. \2A) in MOND/TeVeS. The radius of the 
filament is Rf = 2.5 ft"' Mpc, and the overdensity within the filament is taken as 20 times the mean density po in accordance with the results of Colberg et al. 
(2005). Note that, for consistency, the Newtonian results are based on a flat ACDM cosmology with £2,,, = 0.3 and Q.\ = 0.7. 



In their ACDM large-scale structure simulation, |Colberg 
|et al.'('2005 ) have shown that there are close cluster pairs with 
a separation of 5h^^ Mpc or less which are always connected 
by a filament. At separations between 15 and 20 /;"' Mpc, 
still about a third of cluster pairs are connected by a filament. 
On average, more massive clusters are connected to more fil- 
aments than less massive ones. In addition, the current sim- 
ulation indicates that the most massive clusters form at the 
intersections of the filamentary backbone of large-scale struc- 
ture. For straight filaments, the radial profiles show a fairly 
well-defined radius R j beyond which the profiles closely fol- 
low an r^^ power law, with Rj being around 2.0 hr^ Mpc for 
the majority of filaments. The enclosed overdensity within 
varies from a few times up to 25 times the mean density, in- 
dependent of the filament's length. Along the filaments' axes, 
material is not distributed uniformly. Towards the clusters, 
the density rises, indicating the presence of the cluster infall 
regions. 

As previously stated, we will assume that filamentary struc- 
tures have similar properties in MOND/TeVeS and in a CDM 



dominated universe based on GR. Our assumption is based on 
the /iHDM cosmology (see ^ and on the fact that filaments 
are generic^ and have similar characteristics in hot dark mat- 
ter (HDM) and CDM scenarios ( |Knebe||2003) |Knebe et all 
|2002[[Bode et al. 1200 1*). For instance, neutrinos are known to 
collapse into sheets and filaments in HDM simulations. Con- 
cerning our filament model introduced in Sec. |2] we therefore 
take the filament's radius as /?/ = 2.5 /i"' Mpc, and taking the 
overdensity within the filament to be 20 times the intergalac- 
tic mean density po, we set 5 - 20, with 6 being the density 
contrast defined by 



6^ 



P-Po 
Po 



(22) 



Again, note that choosing the yuHDM cosmology implies that 
filaments do not solely consist of baryonic matter but need 
an additional matter component, i.e. neutrinos, within the 

' Note that the occurrence of filamentary structures is a generic feature of 
gravitational collapse from a Gaussian random field which does not depend 
on the specific form of the law of gravity. 
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Fig. 4. — Density profile p(r) (top left), radial evolution g{r) of the total gravitational acceleration (top right), deflection angle £!'(>') (bottom left) and convergence 
K(y) (bottom right; k = y = (I - A"')/2) in Newtonian (dashed line) and MONDian (solid line) gravity for the oscillating density model given by eq. j26j 
(e = 90°), assuming So = 4, Rf = 2.5 /i"' Mpc, zi = I, Zs = 3, and the flat /iHDM cosmology in eq. j23j for MOND/TeVeS. Note that, for consistency, the 
Newtonian results are based on a flat ACDM cosmology with il,„ = 0.3 and : 



MOND paradigm, which is inferred from the previously men- 
tioned discrepancies between dynamical and visible mass on 
galaxy cluster scale (Sanders 2007) as well as from the need 



for su ch a component to explain the CMB (Skordis et al. 
|2006l l. 

On the other hand, analyzing the Perseus-Pisces segment, 
Milgrom| (| 19971 concluded that a MONDian description of 



filaments would not need any additional non-baryonic mass 
component. Due to rather large systematic uncertainties, how- 
ever, this result remains highly speculative and does not rule 
out our approach where filamentary structures have higher 
densities. Nevertheless, we will also include this case, where 
filaments consist of baryonic matter only, into our analysis. 
Since the absolute density of a filament in this situation is 
approximately by a factor 10 - 100 smaller than in juHDM, 
we do expect the MONDian influence to become even more 



important. Encouraged by the MOND simulation of Knebe 
|& Gibson| (2004), we shall stick to the assumption that both 
shapes and relative densities of filaments are similar to the 
CDM case when considering a universe made out of baryonic 
matter only, thus keeping the choice 6 = 20. 



: 0.7. 

In order to calculate the intergalactic mean density and 
the necessary angular diameter distances for lensing, we still 
need to set up a cosmological model in TeVeS. Depending 
on whether or not we assume massive neutrinos to be present 
in our universe, there will be a different cosmological back- 
ground. If neutrinos are taken into account, then, for sim- 
plicity, we shall use the flat fiHDM cosmology based on the 
parameters of |Skordis et al.| ( [2006| l, 

Q„, = 0.22, Qa = 0.78. (23) 

Considering a universe with baryons only, we choose a flat 
minimal-matter cosmology which is described by 

n,„ = 0.05, Qa = 0.95. (24) 

Furthermore, we shall set h = 0.7 and calculate the model- 
dependent intergalactic mean density po according to 

po ^ n,„pAl + zi)\ (25) 

where pc - 3i/g/87rG is the critical density and zi is the lens 
redshift, i.e. the filament's redshift. Although equation ( |24] i 
has problems in explaining CMB observations due to its pre- 
diction of the last scattering sound horizon, both of the above 
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Fig. 5. — Density profile p(r) (top left), radial evolution g{r) of the total gravitational acceleration (top right), deflection angle £!'(>') (bottom left) and convergence 
K(y) (bottom right; k = y = (I - A"')/2) in Newtonian (dashed line) and MONDian (solid line) gravity for the oscillating density model given by eq. \26\ 
(6 = 90°), assuming So = 4, Rf = 2.5 /i"' Mpc, = 1, = 3, and the flat minimal-matter cosmology in eq. {24) for MOND/TeVeS. Note that, for consistency, 
the Newtonian results are based on a flat ACDM cosmology with Cl,„ = 0.3 and ■ 



cosmological models will be sufficient for assigning the dis- 
tances of lenses and sources at redshifts z < 3 in the context 
of gravitational lensing. However, note that the simplicity of 
these models does not affect the upcoming analysis as we will 
limit ourselves to order-of-magnitude estimates only. 

Concerning the framework of standard Newtonian gravity, 
we shall use a flat ACDM cosmology with Q„, = 0.3 and 
D.A - 0.7, allowing a consistent comparison to our results in 
MOND/TeVeS. 

3.1. The jiHDM Scenario 

Using the cosmological parameters specified in equation 
( [23] l within MOND/TeVeS and considering a filament which 
is inclined by an angle 9 - 90° to the line of sight, both the 
Newtonian and the MONDian deflection angle as well as the 
corresponding convergence are plotted in the bottom left and 
bottom right panel of Figure |2j with the filament placed at 
redshift z/ = 1 and background sources at Zs - 3. Whereas 
the Newtonian signal is rather small, kj^ < 10"^, the fila- 
ment can create a convergence on the order of k ~ 0.01 in 
MOND/TeVeS, even in the outer regions where kn - Q if we 



: 0.7. 

take into account that it can have other orientations, i.e. a 
different inclination angle 0. For example, a nearly end-on fil- 
ament, = 10°, has a lensing power 6 times larger than that 
of a face-on filament, = 90°. 

Using equation ([TO]l, we therefore infer that a single 
MOND/TeVeS filament may generate a shear signal which is 
on the same order as the convergence, y ~ 0.01, as well as 
an amplification bias at a 2% level. A"' ~ 1.02. In addition, 
we present the density p(r) and the radial evolution of the to- 
tal gravitational acceleration g(r) for MOND and Newtonian 
dynamics in the top left and top right panel of Figure |2] re- 
spectively. Note again that, for consistency, the Newtonian 
results are based on a flat ACDM cosmology with Q,„ = 0.3 
and Qa - 0.7. 

3.2. The Baryons-only Scenario 

Now let us switch to the baryonic cosmological background 
given by equation ( p4j i. Keeping all remaining parameters ex- 
actly the same as in the last section, the corresponding re- 
sults are presented in Figure |3] Although the convergence is 
slightly smaller than in the //HDM case, roughly by a factor 
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of 1.5 - 2, we find that also in this case single filamentary 
structures are capable of producing a lensing signal which is 
of the same order, k ~ y ~ 0.01. Again, this is even true 
outside the "edges" of the filament's projected matter density, 
accounting for the fact that the inclination angle may vary, 
0° < 6l< 90°. 

4. OSCILLATING DENSITY MODEL 

Matter density fluctuations are steadily present throughout 
the universe. Thus, as a more realistic approach, we shall 
use a fluctuating density profile to describe a filament and its 
surrounding area including voids, i.e. regions in the universe 
where the local matter density is below the intergalactic mean 
density. To keep our analysis on a simple level, let us write 
the density fluctuation as (r still denotes radial coordinate in 
cylindrical coordinates) 



6o 



0, 









sin — 







r<2R 



(26) 



r>2R 



where 6{r) is the density contrast defined in equation ( [22| i, 
6o - 4is the density fluctuation amplitude (this value ensures 
a positive overall matter density), and Rf = 2.5 /i"' Mpc is 
again the filament's characteristic radius. Multiplying with 
the mean density po and integrating along the radial direction, 
we find that the mass per unit length enclosed by an infinite 
cylinder of radius r reads as (note that we neglect the contri- 
bution due to the mean density background) 



' 2po6oRi 



M{r) 



1 



cos 



R 



r <2R 



(27) 



0, 



r > 2Rf, 

where po is the mean intergalactic matter density given by 
equation p5| ). From equation ( |27] l, we directly see that the 
Newtonian gravitational acceleration in this case is 

gN(r)^——-. (28) 
InL r 

Using equations ([3]), ( [T7] i and ( [28] l, we are now able to nu- 
merically calculate the lensing properties of this configura- 
tion. Choosing lens and source redshift again as z/ = 1 and 
Zs - 3, respectively, and assuming the cosmological back- 
ground models previously introduced in Sec. [3] the result- 
ing deflection angle as well as the convergence are shown in 
the bottom panel of Figures [4] (flat yuHDM cosmology) and [5] 
(flat minimal-matter cosmology), assuming - 90°. Here the 
occurrence of negative /c-values simply reflects the fact that 
our model (|26]l generates a local underdensity, 1 + 6{r) < 1, 
with the overall matter density p being non-negative at any 
radius. Compared to the Newtonian case where < 10""*, 
we again find that a face-on TeVeS filament may cause a sig- 
nificantly larger lensing signal, which is now on the order of 
K ~ y ~ lO""* within both TeVeS cosmologies. As the results 
of the //HDM and the minimal-matter cosmology approxi- 
mately differ by a factor 1.5-2 just as in ^ the order-of- 
magnitude lensing effects caused by TeVeS filaments are also 
in this case more or less cosmologically model-independent. 

Close to the filament's axis, where k ~ A x 10"^, one can 
actually have a lensing signal k - y - 0.01 assuming that 
the inclination angle is small, 9 < 20°. Although such angles 
correspond to rather special configurations, we may conclude 



TABLE 1 

Parameters oe the superimposed EiLAME>iTS in ^ 





PA. 


Incl." 


Shift from Origin' 




Plane 


(deg) 


(deg) 


(kpc) 


Redshift z 


2 


90 


12 


(0,-150) 


0.25 


3 


45 


45 


(600,0) 


0.30 



"Inclination of the filament's axis to the line of sight 

''Shift of the filament's projection in the con'esponding redshift plane 

that also for our simple oscillation model, single TeVeS fila- 
ments can generate a lensing signal ~ 0.01, which is similar 
to our result in ^ However, note that the above discussion is 
based on the choice of equation (|26| and 6q - 4. Considering 
a higher overdensity along its axis, even a face-on filament 
described by a similar fluctuating profile could easily create a 
shear field y ~ 0.01 fory < Rf. 

5. SUPERIMPOSING FILAMENTS WITH OTHER OBJECTS 

To demonstrate the contribution of filamentary structures 
to the lensing map of other objects, e.g. galaxy clusters, we 
superimpose two differently orientated filaments with a toy 
cluster along the line of sight, assuming the previously in- 
troduced yuHDM cosmology and different redshifts for each 
component. If all objects are sufficiently far away from each 
other (> lOOMpc), we may approximately treat them as iso- 
lated lenses at a certain redshift slice, i.e. the corresponding 
deflection angles can be calculated separately**. Thus, we may 
resort to the well-known multiplane len s equation ( ,Blandford] 
& NarayanI 1986} [Schneider et al.|1992| i: 



(29) 



where n is the number of lens planes, dij corresponds to the 
angular diameter distance between the i-th and the j-th plane, 
and ^, is recursively given by 



ii ^^ii-J] 2 < / < «. (30) 



/-I 



Comparing equation ( |29| ) to the lens equation for a single lens 
plane, we identify the total deflection angle, 



(li) = ^i(li) + y ^aiiii) = 4 + Sf, (31) 



where and aj are the deflection angle of an isolated cluster 
at zi and an additional contribution due to the superimposed 
filaments, respectively. Analog to the case of a single plane, 
further lensing quantities such as the total convergence and 
the total shear can be calculated from equation ( [3T| ), using the 
general relations introduced in ^ For simplicity, we shall 
assume that the clust er's TeVeS potential follows the "quasi- 
isothermal" profile of jAngus et al.| ( 2007) l: 



(D(/0 = v2log Jl + 



\r - ro\- 



(32) 



where v is the asymptotic circular velocity, /? is a scale length, 
and ro is the center's position. 

^ Note that, in general, one would have to solve the full non-linear TeVeS 
scalar field equation, which is beyond the scope of this work. 
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Fig. 6. — Superposition of two filaments with a toy cluster along the line of 
structures along the line of sight (top left) and the filaments' contribution Ak = 
reduced shear, Agi = 7m,,i/(1 - k,o,) - yc.i/(i - Kc) and Ag2 = ytot.lKX - "wi) 

Concerning the numerical setup, we set = 2 x lO*" 
(km s"')^ and p - 200 kpc, fixing the cluster's redshift to 
Zi - 0.2. Furthermore, we choose the uniform filament model 
discussed in ^ and assume that filaments have a constant 
overdensity of 5 = 20 as well as the same characteristic ra- 
dius Rf = 2.5 /z"' Mpc. While the cluster is centered at the 
origin - - 0), the two filaments are set up according to 
the parameters given in Table [T] Finally, we place the source 
plane at a redshift of Zs = 1 ■ Note that this specific setting cor- 
responds to a more realistic lensing configuration compared to 
our analysis in the sections above, with our choice again being 
motivated by results based on a ACDM universe. 

From the top right panel of Figure [6] we see that the fil- 
aments' contribution to the total convergence map, Ak - 
Kior - Kc {kc is the cluster's convergence map in absence of any 
filamentary structures along the line of sight) is comparable to 
our previous findings, with the signal again being on the order 
of 0.01. Also, note the distortion effects caused by the clus- 
ter and the peak close to the region where the two filaments 
overlap. Obviously, the contribution pattern depends on the 
actual configuration as well as on the type and amount of the 
considered objects along the line of sight, and can generally 
be quite complex. In addition, we present the changes in the 
reduced shear components, Agi - 'yrot.\/{^-'<rot)-7c.i li^-Xc) 
and Ag2 - ytot,2Kl - Xtor) - rf,2/(l - Hc), due to the filaments' 
presence in the bottom panel of Figure |6] 

At this point, we should emphasize that we have consid- 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

[Mpc] 
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sight; shown are the cluster's convergence map in absence of any filamentary 
Ktot - Kc to the total convergence {top right), as well as to the components of the 
~ 'yc,2/(l - Kc), respectively (bottom panel). 

ered the impact of filamentary structures alone. Depending 
on their particular position along the line of sight, additional 
(foreground) objects such as galaxies, galaxy clusters and/or 
voids might (locally) contribute on a comparable level or even 
exceed the signal caused by filaments. Of course, this fur- 
ther complicates the interpretation of the corresponding lens 
mapping, and we conclude that, in general, extracting the fil- 
aments' contribution can pose quite a challenge. 

6. CONCLUSIONS 

In this work, we have analyzed the gravitational lensing ef- 
fect by filamentary structures in TeVeS, a relativistic formula- 
tion of the MOND paradigm. 

For this purpose, we have set up two different cosmolog- 
ical models in TeVeS: the so-called /iHDM cosmology in- 
cluding massive neutrinos on the order of 2 eV which have 
already been proposed as a remedy for the discrepanc ies be- 
tween dynamical and visible mass on cluster scales (Sanders 
[2007l l as well as for the CMB ( |Skordis et a l. 2006), and a 
simple minimal-matter cosmology accounting for a universe 
which is made up of baryons alone. Encouraged by several 
HDM simulations and the fact that filamentary structures are 
generic, we have assumed that the properties of such struc- 
tures, i.e. their shape and relative densities, are similar in 
CDM and MOND/TeVeS scenarios independent of the par- 
ticularly used cosmological background. 

Modeling these filaments as infinite uniform mass cylin- 
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ders, we have derived analytic expressions for their lensing 
properties in MOND/TeVeS and Newtonian/GR gravity. Re- 
gardless of the actual used cosmological background, we have 
shown that TeVeS filaments can account for quite a substantial 
contribution to the weak lensing convergence and shear field, 
K ~ J ~ 0.01, as well as to the amplification bias. A"' ~ 1.02, 
which is even true outside but close (y ~ 2R j) to the projected 
"edges" of the filament's matter density. Exploring a sim- 
ple oscillating density model of a filament and its surround- 
ing area, we have found that the lensing signal in this case is 
generally smaller, but can still be of the same order, taking 
into account that the filamentary structures may be inclined 
to the line of sight by rather small angles (6 < 20°). In addi- 
tion, we find that there is fundamental difference between GR 
and MOND/TeVeS considering idealized cylindrically sym- 
metric lens geometries: wherever the projected matter den- 
sity is zero, there will be no distortion as well as no ampli- 
fication effects, i.e. image and source will look exactly the 
same. In the context of MOND/TeVeS, however, this changes 
as one can have such effects in these regions. Finally, we have 
demonstrated the impact of filaments onto the convergence 
map of other objects by considering superposition with a toy 
cluster along the line of sight. Again, our results have shown 
an additional contribution comparable to that of a single iso- 
lated filament and that the contribution pattern of filaments 
can be generally quite complex. 

Here we have considered the lensing signal generated by 
single filaments alone. Simulating the cosmic web in a stan- 



dard ACDM cosmology, Dol ag et al.| ( 2006 1 have found a 
shear signal y ~ 0.01 - 0.02 along filamentary structures, 
which seems quite similar to what MOND/TeVeS can do. 
Note, however, that this signal is entirely dominated by the 
simulation's galaxy clusters, with the filament's signal being 
much smaller, approximately on the order of 10""* - lO^^*. 

Although our analysis is mainly of theoretical interest, the 
above result points to an interesting possibility concerning 
recent measurements of weak lensing shear maps. For in- 
stance, the weak shear signal in the "dark matter peak" of 
Abell 520 ([Mahdavi et al.||2007|l is roughly at a level of 



0.02, which is comparable to what filaments could pr oduce in 
MOND/TeV eS, but not in Newtonian gravity (also cf. Famaey 
et al. ( 2008| ). Therefore, we conclude that filamentary struc- 
tures might actually be able to cause such anomalous lensing 
signals within the framework of MOND/TeVeS. 

In principle, the predicted difference in the weak lensing 
signal could also be used to test the validity of modified grav- 
ity. As several attempts to detect filaments by means of weak 
lensing methods have failed so far, e.g. the analysis of Abell 
220 and 223 by Dietrich et al. ( 2005 1, this might already be 
a first hint to possible problems within MOND/TeVeS grav- 
ity. On the other hand, shear signals around y ~ 0.01 are still 
rather small to be certainly detected by today's weak lensing 
observations, and lacking A^-body structure formation simu- 
lations in the framework of MOND/TeVeS, we cannot even 
be sure about how filaments form and how they look like in 
a MONDian universe compared to the CDM case. Clearly, 
more investigation is needed to gain a better understanding 
about the impact of filamentary structures. 
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